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Abstract 

Large-scale cosmic microwave background anisotropies in homogeneous, glob- 
ally anisotropic cosmologies are investigated. We perform a statistical analysis 
in which the four-year data from the Cosmic Background Explorer satellite 
is searched for the specific anisotropy patterns predicted by these models 
and thereby set definitive upper limits on the amount of shear, {a/H)^ and 
vorticity, {lo/H)q, which are orders of magnitude stronger than previous con- 
straints. We comment on how these results might impact our understanding 

of primordial global anisotropy. 
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Typeset using REVT^ 



Fluctuations in the cosmic microwave background radiation (CMBR) provide us with 
a clean and unique probe of the structure of our universe on both small and large scales 
||1|] . Current experiments are allowing us to learn about processes in the early universe with 
unparalleled precision, and the prospect of future terrestrial and satellite missions give us 
reason to hope for a clear and accurate picture of the universe in the near future. The main 
working hypothesis in our attempts to understand the universe is that we live in what was 
originally a homogeneous and isotropic spacetime, the Friedmann-Robertson- Walker (FRW) 
cosmology. The current belief is that some physical mechanism generated perturbations (ei- 
ther primordial, with inflation, or actively, with defects) that evolved through gravitational 
collapse to form the structures we see now. The smoothness of the CMBR seems to be 
consistent with such a picture. 

Inflation provides the primary motivation for believing that, at least in an initially homo- 
geneous space-time, any primordial anisotropy has been removed. However, inflation is by 
no means a generic initial condition for the universe. In the absence of inflation, one is likely 
to have commenced with global metric perturbations that while eventually seeding structure 
formation may result in a universe that is only asymptotically Friedmann-Robertson- Walker. 
Bianchi models provide a generic description of homogeneous anisotropic cosmologies. In the 
spirit of studying an alternative to inflation, we study below the experimental constraints 
if our universe is only asymptotically FRW. There are distinct features depending on the 
overall geometry and homogeneity class of the model , and in a pioneering paper, Collins 
and Hawking used analytical arguments to find upper bounds on the amount of shear (ctq) 
and vorticity (ujq) in the universe today, from the absence of any detected CMBR anisotropy. 
A detailed numerical analysis of such models |^ used experimental limits on the dipole and 
quadrupole to refine limits on universal rotation. More recently. Barrow has argued that 
there is no "isotropy problem" in such cosmological models which are maximally anisotropic 
(i.e. in which shear and vorticity have decayed only logarithmically since the Planck time): 
if one assumes an equipartition of energy between all modes of the gravitational interaction 
at the Planck time then the present amplitude of CMBR fluctuations should be compatible 
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with current observational limits. 

In this Letter we improve on previous bounds on the total shear in the universe. We 
use the most recent data from the Differential Microwave Radiometers aboard the Cosmic 
Background Explorer (COBE) ^ to constrain the allowed parameters of a Bianchi model of 
type Vllh] this model is asymptotically close to an open FRW universe and has the richest 
anisotropy structure of the models we could consider. As pointed out by Barrow, it is also 
an example of a homogeneous cosmology where the decay in a from the Planck time to the 
present is minimal. 

One can describe Bianchi cosmologies in terms of the metric 

gi„y = -n^n^ + a^[exp(2/?)]ABe;^ef , (1) 

where Ua is the normal to spatial hypersurfaces of homogeneity, a is the conformal scale 
factor, Pab is a 3 matrix only dependent on cosmic time, t, and are invariant covector 
fields on the surfaces of homogeneity, which obey the commutation relations 

The structure constants C^^ can be used to classify the different models. We shall focus on 
type VII/i which has structure constants C^^ = C21 = 1, C21 = C^^ = \/h. It is convenient 
to define the parameter x = \Jh/l — fio, which determines the scale on which the principal 
axes of shear and rotation change orientation. By taking combinations of limits of VL and x 
one can obtain Bianchi I, V and VIIq cosmologies. 

We are interested in large-scale anisotropics so it suffices to evaluate the peculiar redshift 
a photon will feel from the epoch of last scattering [Is) until now (0) 

AT^(f) = {eui), - {eu,)u - f'f^f'ajkdT (3) 

J Is 

where f = (cos 6 sin (f), sin 6 sin 0, cos (p) is the direction vector of the incoming null geodesic, 
u is the spatial part of the fluid four- velocity vector and to first order, the shear is aij = drPij. 
To evaluate expression (j^), one must first of all determine a parameterization of geodesies 
on this spacetime. This is given by 
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tan(^^) = tan(y) exp[-(r - To)Vh] 
9{T) = 9o + {r-ro) 



^ \n{sm\^) + cos\^) exp[2(r - To)Vh]} (4) 



Solving Einstein's equations (and assuming that matter is a pressureless fluid) one can 
determine u and 0"^. A general expression for was determined in [Q: 



X { [sin 00 COS - sin COS 6';s ( 1 + ) ] 

- IZ ^ 20[cos(^) + sm(^)]-j-L_} 

(5) 

As shown in the "patterns" in such a model are easy to describe: for i^o < 1 and 

a finite x one will obtain a spiral with approximately = 2/7tx complete twists, focused 
towards the axis of rotation with an angular size of order I/Qq. Taking x oo will leave 
only a hotspot. 

We will approach the problem of constraining the parameters of these Bianchi models 
in the following way. For fixed values of the parameters x and Qq, we will attempt to place 
upper limits on the amplitude of the shear {a/H)o (or equivalently vorticity {uj/H)o 0). 
The statistics problem we face differs substantially from the situation encountered in placing 
constraints on more standard models. In standard cosmological models, the predicted CMB 
anisotropy is a realization of an isotropic Gaussian random field, and its statistical properties 
are therefore entirely characterized by the power spectrum C/. In contrast, in the Bianchi 
models at least part of the CMB fiuctuation comes from the large-scale anisotropy of space; 
this contribution to the anisotropy takes the form of a particular pattern on the sky, and 
is not described by the statistics of a Gaussian random field. We therefore require different 
statistical techniques from those used in previous analyses. 

Let AT(r) be the temperature fluctuation in the direction of the unit vector r. We 
assume that AT is the sum of two contributions: 
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Ar(f) = AT^(r) + AT,(r). (6) 

Here AT4 is defined in (||) and AT/ represents the "isotropic" residual fluctuation caused 
by variations in the density and gravitational potential. We call AT/ "isotropic" because we 
assume that it is described by the statistics of an isotropic Gaussian random fleld. That is, 
if we expand it in spherical harmonics, 

AT,(f) = 5]a,„l^„(f), (7) 

l,m 

then the coefficients aim are independent Gaussian random variables with zero mean and 
variances given by an angular power spectrum Ci: {aimdi'm') = Ci^ii'^mm'- 

One can naively assume that this set of Gaussian perturbations was generated through 
the ampliflcation of quantum fluctuations, as in the case of inflation (however see 0) and it 
is quite conceivable that this initial set of perturbations is strongly biased (or anti-biased) 
due to the large anisotropy at early times. To change our results considerably, this would 
have to compensate (anticompensate) late time evolution of the overall anisotropy on many 
scales, i.e. there would have to be a strong correlation between the primordial quantum 
generation of perturbations and subsequent classical evolution of the different temperature 
variables such that they would interfere (destructively or constructively) for many modes of 
the temperature autocorrelation function, something that we believe to be unlikely. Another 
possibility is that subsequent evolution of perturbations will be locked in to the speciflc 
orientation of the large scale anisotropy, but if we assume that the dominant source of 
perturbations are scalar and therefore are only sensitive to the overall volume change of 
the spacetime, than we can discard this hypothesis. This reasoning does not hold if the 
dominant source of perturbations comes from tensor modes, but we shall not consider this 
possibility here. 

The anisotropic component AT^ does not obey Gaussian statistics. Rather, for flxed 
values of x, Qq, and a, we can compute the exact pattern A{r) of the CMB anisotropy. 
Unfortunately, we do not know the orientation of this pattern on the sky. We can say that 

AT^ = A(Rr), (8) 
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where A is the known pattern of anisotropy and R G 5*0(3) is a 3 x 3 rotation matrix. The 
matrix R can be specified by three Euler angles {6, ijj, (p), but we have no knowledge a priori 
of the values of these angles. 

We will place constraints on the shear following standard frequentist statistical practice. 
(This is in contrast to the Bayesian philosophy adopted in much cosmological data analysis.) 
For any particular model, we define some goodness-of-fit statistic rj that depends on the data 
(To give a simple, familiar example, when one is trying to estimate the mean of a set of 
data, it is customary to choose rj to be the chi-squared of the data). Having chosen our 
goodness-of-fit statistic, we compute its value r^* using the actual data. We then compute 
the probability P(?7 < 77^,) that a random data set would have produced a value as good a 
value as the actual data or better. If this probability is large, then we say that the model is 
inconsistent with the data. It is customary to choose a significance level Pq, say 0.95, and 
say that a particular model is ruled out at that significance level if P{ri < rj^) > Pq. 

Our choice of rj is as follows. Each pixel di of our data set contains contributions from 
both intrinsic CMB anisotropy and noise: 



here B represents the DMR beam pattern P], rj is the direction on the sky of the ith data 
point, and the star denotes a convolution. According to our model, AT includes the two 
contributions shown in equation (|^). Ni is the noise in pixel i. For the COBE data, it is 
an excellent approximation to take to be Gaussian with zero mean. The correlations 
between the noise in different pixels are negligible 0, so {NiNj) = crfdij. 

Before we give the actual definition of the statistic 77, let us consider a simpler case. 
Suppose that we knew the geometrical parameters a;, VLq, and a, and the rotation matrix 
R that defines the orientation of the pattern A on the sky. Then the anisotropic part AT4 
of the CMB fiuctuation would be completely specified, but the isotropic portion AT/ and 
the noise Ni would be completely unknown. In this situation, we could define a natural 
goodness-of-fit statistic in the following way. Let Aq be the noise-weighted mean-square 




(9) 
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value of the data: 



= E I- (10) 

Now let Af be the mean-square value of the residuals after we have subtracted off the known 
anisotropic portion: 

^,_^ id,-iATA^B)ir,)f ^ ^^^^ 

If our model is correct, then we expect to be smaller than Aq: if we have correctly 
removed a portion of the signal, then the residuals should be smaller, on average. On the 
other hand, if our model is incorrect, then attempting to remove the anisotropic portion 
should increase the residuals. The difference between Aq and A^ is therefore a natural 
choice of goodness-of-fit statistic. In practice, it is more convenient to divide by Aq, in order 
to make the results more weakly dependent on the amplitude of the isotropic cosmological 
signal AT/. We therefore define our statistic to be 

m = ^2-^. (12) 

^0 

In fact, of course, we cannot use the statistic (|T2]), because we do not know the parameters 
necessary to determine ATa- In particular, we do not know R or cr. (We have chosen to set 
up the problem as one of constraining a for fixed x and Qq, so we can assume that we know 
the latter two parameters.) But we can define a new statistic r] whose value is the minimum 
of ?7i over the unknown parameters: 

1] = min ?7i. (13) 

The statistical task we have set ourselves is simple in principle, although it is somewhat 
cumbersome computationally. For fixed values of the parameters x and Qq, we must compute 
the value rj^ of the statistic (|1^) for the real data. We must then perform Monte Carlo 
simulations to determine the theoretical probability distribution of rj to see how consistent 
the actual value is with each theoretical model. We must perform these simulations for a 



variety of different values of the shear a in order to see which values of a are consistent with 
the data. 

Each calculation of r] involves a minimization in a four-dimensional parameter space. 
Since we need to compute r] repeatedly in our simulations, it is important to perform this 
calculation efficiently. We chose to reduce the numbers of pixels in the COBE data set by 
binning pixels together in groups of four {i.e., working in "pixelization level 5" rather than 
level 6). Since the anisotropy pattern /ST a tends to have power on larger scales than either 
the noise or the isotropic signal, this binning does not reduce our sensitivity very much. 

The task of finding the minimum in equation ( |13|) is not trivial, since the function 



T]i{a,6,ip,(f) has numerous local minima. We use Powell's method for finding minima, but 
we have to try multiple starting points in order to be confident that we had found, if not 
the true minimum, at least a local minimum that was almost as low as the true minimum. 
We chose to adopt the following procedure. We choose p random points in parameter space 
and evaluate rji at each. Starting from the point that gave the lowest value of rji, we use 
Powell's method to find a local minimum. We then repeat this entire procedure q times, and 
we take the lowest value found to be our statistic 77. After some experimentation, we found 
that choosing p = 10 and g = 4 gave reasonable results. Of course, it is essential to use 
precisely the same procedure for determining 77 in both the real data and the simulations. 

We perform our analysis on the four- year DMR data set IQ. We use the ecliptic pix- 
elization of the data. Before performing any analysis, we average together the two 53 GHz 
maps and the two 90 GHz maps to make a single sky map. The averaging is performed with 
weights inversely proportional to the squares of the noise levels, in order to minimize the 
noise in the average map. In order to reduce Galactic contamination, we excise all pixels 
that lie within the "custom cut" described by the COBE group [|T^ ; this reduces the number 
of pixels in the map from 6144 to 3890. We then remove a best-fit monopole and dipole 
from the map. As mentioned above, we degrade the map from pixelization level 6 to level 
5, and compute the statistic 77 for a grid of points in the Qq-x plane. 

In order to convert these t] values into constraints on the shear, we need to determine 
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the probability distribution of t] via Monte Carlo simulations. We performed simulations on 
a sample of 10 models in the Qq-x plane, using three or four values of a for each model. For 
each choice of the three parameters [Qq, x, a), we created between 200 and 500 random DMR 
data sets according to equation (|^). For simplicity, we assumed that the isotropic component 
AT/ of the anisotropy was given by a scale-invariant power spectrum Cf^ oc /(Z+l), although 
our final results are not very sensitive to this assumption. [Specifically, if we steepen the 
power spectrum to an effective n = 1.5 spectrum (see, e.g., ||l[| for a definition), the limits in 
Figure 1 change by ~ 20%.] We processed each sky map in the same way as the real data 
to determine a value of 1]. 

We found that in every case the probability distribution of rj was slightly skew-positive 
and had tails that were consistent with exponential distributions. We determined the first 
three moments of each probability distribution, and found that each distribution was very 
well approximated by a stretched, offset chi-squared distribution, where the three parameters 
of the distribution (stretch, offset, and number of degrees of freedom) were chosen to fit the 
three moments ||rT|. For points in parameter space where we have not performed simulations. 



we assume that the probability distribution is also well approximated by a stretched chi- 
squared distribution, and we determine the three parameters of the distribution by smoothly 
interpolating between the points where we have performed simulations. 

Having estimated the probability distribution of rj for the various theoretical models in 
this way, we are able to set limits on the shear. For each point in our grid in the Qq-x plane, 
we determine the range of values of a such that P{r] < 7]^) < 0.95. We find that a = is 
always allowed at the 95% confidence level; i.e., we do not detect shear at this level. Figure 
1 shows the upper limits we can set on [a / H)q and {uj/H)q as a function of VLq and x: for 
1^0 = 1 universes we find that {a / H)q < 3 x 10~^ (or {uj/H)q < 10^^) for any x > .05, while 
for fio < 1 the upper bounds are even tighter. 

These values are to be compared with the constraints from [Q] which are typically one 
to two orders of magnitude higher and relied entirely on the quadrupole: at that time 
g ~ 7 X IQ-^K compared to Q ^ 1 - 2 x 10~^K from the COBE DMR data Moreover, 
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FIG. 1. Upper bounds on shear and vorticity: contours of equal {(t/H)q x 10^ and logxo(w/-ff)o 
are shown in the upper and lower panels respectively for a class of Bianchi VII/i models. 

in discarding information from higher moments, their comparison was not sensitive to the 
small-scale structure present in anisotropic models that is associated either with the spiral 
pattern (which introduces power on smaller scales as x increases) or with geometrical focusing 
when i7 < 1. 

Our tighter constraints rule out the Planck equipartition principle for primordial global 
shear and vorticity in its most general form [^]. If we consider logarithmic decay of shear 
through the radiation era due to collisionless stresses, then a rough estimate gives {a/ H)pi ~ 
(cr/if)o(l + Zcurv){,^ + Zf,q)~^{l + \n.{teq / tpi)}^^ whcrc Zcurv (^eg) IS the rcdshift of curvature 
(matter) domination. For the "best" case of i^o = 1 "we obtain {a/H)pi ~ 10^^ — 10^'*. Our 
argument applies to the most general allowed set of globally anisotropic models: generic 
open or fiat homogeneous, anisotropic but asymptotically Friedmann models. While it is 
possible that for more rapid decay of shear one can allow some relic anisotropy, we conclude 
that any residual shear or vorticity is constrained to be so small that it is most unlikely to 
have had a Planck era origin. 
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